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Abstract 

Let {Rt, 0<i<l}bea symmetric a-stable Riemann-Liouville process with Hurst 
parameter H > 0. Consider a translation invariant, /3-self-similar, and p-pseudo-additive 
functional semi- norm ||.||. We show that if H > (3 + 1/p and 7 = (H — (3 — 1/p)^ 1 , then 

lime 7 log P[||i?|| <e] = -K G [-00, 0), 

with /C finite in the Gaussian case a = 2. If a < 2, we prove that /C is finite when R is 
continuous and H > f3 + 1/p + I /a. We also show that under the above assumptions, 

lim£ 7 logP[||X|| < e] = —JC G (-00, 0), 

ej.0 

where X is the linear a-stable fractional motion with Hurst parameter H G (0,1) (if 
a = 2, then X is the classical fractional Brownian motion). These general results cover 
many cases previously studied in the literature, and also prove the existence of new 
small deviation constants, both in Gaussian and Non-Gaussian frameworks. 

Resume 

Soit {Rt, < t < 1} un processus de Riemann-Liouville a-stable symctrique avec 
parametre de Hurst H > 0. Considerons une semi-norme fonctionncllc ||.|| invariante 
par translation, /3-autosimilaire et p-pseudo-additive. Nous montrons que si H > /3+l/p 
et 7 = (H- (3 - l/p)- 1 alors 

lime 7 logP[||i?|| <e] = -K G [-00, 0), 

ej.0 

avec K, finie dans le cas gaussien a = 2. Lorsque a < 2, nous montrons que K est finie 
quand R est continu et H > (3+1 /p+ 1/a. Nous montrons aussi que sous ces hypotheses 

lime 7 logP[||X|| < e] = —K G (-00, 0), 

eiO 

ou X est le mouvement fractionnaire lineaire a-stable avec parametre de Hurst H G (0, 1) 
(lorsque a — 2, X est le mouvement brownien fractionnaire usuel). Ces resultats 
generaux recouvrent de nombrcux cas precedemment etudies dans la litterature et prou- 
vent l'existence de nouvelles constantes de petites deviations, aussi bien dans le cadre 
gaussien que non gaussien. 

Keywords: Fractional Brownian motion - Gaussian process - Linear fractional stable mo- 
tion - Riemann-Liouville process - Small ball constants - Small ball probabilities - Small 
deviations - Stable process - Wavelets. 

MSC 2000: 60E07, 60G15, 60G17, 60G18 
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1 Introduction 



Let X be a random process whose sample paths belong to some functional normed space 
(J-, || ■ ||). The investigation of the small deviations (or small ball probabilities) of X deals 
with the asymptotics of 

W[\\X\\<e] 

when e [ 0, and has proved to be a difficult problem with increasing number of applications 
in Probability, Analysis, Complexity ... etc. We refer to the recent surveys |25| |26j for a 
detailed account on this subject. 

In the literature, this problem is usually studied for a particular class of processes and 
under a particular norm. It remains a great challenge, a kind of "mission impossible", to 
find some principle describing small deviations for general classes of processes and norms, 
rather than tackle the problem case by case. 

The unique successful attempt in this direction was made by W. Stolz [10], who 
obtained estimates 

-oo < liminfe 7 logP[||X|| < e] < limsupe 7 logP [\\X\\ < e] < 0, (1) 

where {Xt, 0<i<l}isa Brownian motion (or more generally a continuous Gaussian 
process with covariance function similar to that of fractional Brownian motion), and 7 
a finite positive parameter depending on the behavior of the (semi-)norm || • || on linear 
combinations of Schauder functions. Many classical semi-norms fell into Stolz's scope: L p - 
norms, Holder and Sobolev semi-norms, Besov norms... etc, and his estimates resumed from 
a general point of view a lot of situations which were previously studied (see the references 
quoted in [39] |30]'). 

The next important issue is the existence of the limit in (^Q), and this is the matter of 
the present paper. Our main result says that if ||.| is a translation invariant, /3-self-similar 
and p-superadditive functional semi-norm - see Definition 1 for more precisions about these 
notions, and {Rt, < t < 1} is a symmetric ce-stable Riemann-Liouville process with Hurst 
parameter H > (3 + 1/p (R can be viewed as a fractionally integrated symmetric a-stable 
Levy process, see Section 2 for a precise definition), then 

lime 7 log P[||.R|| <e] = —K, £ [-00, 0) (2) 

with the rate 7 = (H — (3 — l/p)~ l . From the technical point of view, the main ingredient 
of the proof is a stochastic superadditive inequality which is based upon the extrapolation- 
homogeneity of R, and is then easy to combine with an exponential Tauberian theorem and 
the standard subadditivity arguments. 

Our framework has two secondary, but non-negligible advantages with respect to [4Uj . 
Non-Gaussian stable processes are included, as well as smooth Gaussian ones like integrated 
Brownian motion of arbitrary parameter. 

A special effort is needed to prove that the constant /C in (J2J) is finite. In order to get the 
suitable lower bound on the small deviation probabilities, we use roughly the same method 
as in [22] [101; except that we decompose R along Daubechies' wavelet bases, since Schauder 
system not smooth enough when H becomes too large. 

In the last part of this paper, we extend the above results to a class of self-similar processes 
with long range dependence, the so-called unilateral linear fractional stable motions. These 
processes can be viewed as a possible generalization of fractional Brownian motion with an 
underlying stable noise. Thanks to a nice argument essentially due to of Li and Linde [22] . 
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the problem is reduced to a study of the Schauder decomposition of the associated long 
memory process. 

We conclude this article with a brief survey of concrete results. It seems that our theorems 
synthesize everything that is known about the existence of finite small deviation constants for 
continuous fractional processes under translation-invariant semi-norms. Last but not least, 
we can harvest a bunch of new constants, both in Gaussian and Non-Gaussian situations. 

Still, a major drawback of the above wavelet methods is that they exclude discontinuous 
processes. As a rule, proving lower bound probabilities for processes with jumps requires 
completely different discretization techniques [HJ @2]) and we have no idea how much time 
it will take to enclose them into a global strategy. 

2 Preliminaries 

2.1 Parametrization: (/9,p)-semi-norms 

Let X be the set of all closed bounded intervals of R. Consider T a linear space of functions 
from R to R and, for each I El, let Ti be a linear space of functions from I to R such that 
fj E Tj for every / G JF, where fj stands for the restriction of / to /. 

We define a semi-norm || • || on T as a family {|| • ||j, I G X} of functionals mapping T\ to 
R+ such that 1 1 A/| |, = |A|| |/| |, and ||/ + <?|| 7 < \\f\\j + \\g\lj for every A G R,f,g E ft. We 
will use the notation \\fWj = \\fi\\j for every / G T, I G X. In the remainder of this paper 
we will assume that || • || satisfies the following assumptions, which are verified by all the 
classical semi-norms: 

(A) || • ||j < || • \ \j for every I, J El such that I C J. (Contractivity) 

(B) ||/||j_ c = ||/(- — c)\\j for every / G J 7 , I G X and c G R. (Translation-invariance) 

Definition 1 Let j3 G R, p G (0, +oo] and \ \ ■ \\ be a contractive and translation-invariant 
semi-norm on T . We say that \ \ ■ \ \ is an upper p)-semi-norm if it satisfies the following 
properties: 

(C) ||/(c-)||j/ c = c 13 H/llj for every f E T, I E X and c> 0. (/?-self-similarity) 

(D) For every ao < . . . < a n G R and f G T 

H/ll[ao,a„] ^ U l/l '[«o.«x] + "" + ll / ll[an-i,«n]J ^ P < +°°, 

< (p-superadditivity) 

, ll/ll[ao,an] ^ SU p(ll/ll[ao,a 1 ]'---H/llK_ 1 ,a n ]) if V = +°°. 

In the following we will denote by U(/3,p) the set of upper (f3,p)-semi- norms, and set U 
for the union of all U(/3,p)'s. 

Remarks (a) Of course it suffices to take n = 2 in the definition of p-superadditivity. We 
wrote the property in this form in order to make it symmetric with the p-subadditivity and 
the corresponding lower (/3,p)-semi-norms, which will appear just below. 

(b) From the inequality (a + b) q > a q + b q for every a,b > and q > 1, it follows that 
U(P,p) C U(/3, p') whenever p' > p. 

(c) In the definition of upper p)-semi-norm, one can ask for the possible values of the 
parameters f3 and p. By contractivity, homogeneity and translation-invariance, it is easy 
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to see that if < ||l||[o,i] < co (where 1 stands for the unit function), then necessarily 
(5 + 1/p < 0. However, the examples (c)-(f) below show that this inequality is no more true 
whenever ||l||[o,i] = 0. 

We stress that most of the usual semi- norms belong to U: 

Examples (a) The supremum semi-norm, which is given by 

11/11/ = sup 
tei 

for every / el, belongs to U(0, +oo). 

(b) The Lp-semi-norm, p > 1, which is given by 

i/p 



xfi \f(s)\ p ds 

for every I El, belongs to U(— l/p,p). 

(c) The //-Holder semi- norm, < rj < 1, which is given by 

J/(*)-/(*)l 

s<tei 



i = sup 



for every I El, belongs to U(r/,+oo). In particular, the oscillation semi-norm {rj = 0) 
belongs to U(0, +oo). Similarly, the Calderon-Zygmund semi-norm 

ii/ii, = SUP \m±m^m 

3 <tei it ~ s ) 

belongs to U(l,+oo), and the 77-Lipschitz semi-norm, rj > 1, which is given by 

11/1,,= BUP '^w-ff wl 

s<tei \t — s\' 

where n < rj < n + 1 (and by its Calderon-Zygmund analogue for rj = n), belongs to 
Ufa, +oo). 

(d) The (strong) p- variation semi-norm, p > 1, which is given by 

/ n \ VP 

ii/ii/ = sup £wwfe-i)n 

\to<-<*r.e/ i=1 / 
for every I El, belongs to \J(0,p). 

(e) The (?y,p)-Sobolev semi-norm, p > 1 and < rj + 1/p < 1, which is given by 

p \ i/p 



\f(t)-f(s)\ 

1 " U/7/l li-si^vp 



(is (ft 



for every I El, belongs to U(r/ — l/p,p). 

(f) The (t/, p, g)-Besov semi-norm, 77 > and p, q > 1, which is given by 
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for every I el, where 



u P) i{t, f) = sup I / \f(x - h) - f(x)\ p dx 

\h\<t \JI h 

Ih = {x G / 1 x — h G 1} and |/| stands for the Lebesgue measure of /, belongs to \J(rj — 
1/p, +00). Notice that the usual (i],p, q)-Besov semi-norm, where the L p -semi-norm is added 
in the definition (see e.g. JU), does not belong to U since L p -semi- norms have a different 
self-similarity index. 

It follows from Remark (b) above that the concept of upper ((3, p)-semi-norm is not sharp 
enough for our further purposes. We need some kind of inverse property to (D), in order to 
work with a family of disjoint sets of semi- norms. Let C l K (resp. C^) denote the class of all 
/-times continuously differentiable (resp. continuous) functions with compact support. For 
technical reasons, we may have to make the following assumptions on 1 1 • 1 1 which, again, are 
verified by all the classical semi-norms: 

(E) 3 / G N such that C l K C T. (Smooth-finiteness) 

(F) For every I el and /, f n G Ti, 

fn~*f uniformly on I ||/||/ < limsup ||/ n ||/- (Lower semi-continuity) 

Remark Clearly, /-smooth-finiteness yields ||/||/ < +00 for every / of class C l and every 

I el. 

Definition 2 Let (3 G R, p G (0, +00] and \\ ■ \\ be a contractive, translation-invariant, 
smooth- finite, lower semi- continuous and (3 -self-similar semi-norm on J-. We say that \ \ ■ \ \ 
is a lower (/3,p)-semi-norm if it satisfies the following property: 

(G) There exists C p > such that for every ao < ... < a n G R and f G C l K verifying 
/(a ) = • • • = f(a n ) = 0, 

' H/IUan] ^ ^(11/1^] + ••• + H/llfa„_ 1 , a „]) 1/P if P <+00, 
* (p-subadditivity) 

, \\f\\ [ao ,a n] < Ox, SUP (il/Iltoo^], - - - . Il/ll[«_x,-]) if P=+~- 

In the following we will denote by L(/3,p) the set of lower (/3,p)-semi-norms, and set L 
for the union of the L(/3,p)'s. We will say that || ■ || is a {(3, p)- semi-norm if it belongs to 
N(/3,p) = XJ(P,p) n li(f3,p) and analogously we will set N for the union of N(/3,p)'s. 

Remarks (a) Since the constant C p may be larger than 1, it is not enough to take n = 2 
in the definition of p-subadditivity, contrary to p-superadditivity. 

(b) From the inequality (ao + • • • + a n ) q < + • • • + a^, for every ao, . . . ,a n > and 
< q < 1, it follows that L(/3,p) C L(/3,p') whenever p' < p. 




5 



(c) Again, one can ask for the possible values of the parameters [5 and p in the definition of a 
lower (p, /3)-semi-norm. By contractivity, homogeneity and translation-invariance, it is easy 
to see that when the inequalities in (G) hold for the function / = 1 and for all partitions, 
then necessarily (3+ 1/p > 0. This inequality is actually true for all available examples, but 
we got stuck in proving that in full generality. 

(d) It is not difficult to see that N(/?,p) n N(/3',p') = as soon as (/3,p) ^ (f3',p'). 

Notice that each of the above examples (a)-(d) belongs to L, and hence to N. Actually 
in each case we had chosen the smallest possible parameter p, i.e. we could have written 
|| • || £ L(/3,p) as well as || • || £ TJ((3,p). This fact is trivial for examples (a)-(c) where we 
can take C p = 1, a bit more involved in example (d) where we have to take C p = 2 1_1//p as 
soon as p > 1 (notice that here the condition /(ao) = • • • = f(a n ) = is essential). 

However, (n,p, q)-Besov and (r/,p)-Sobolev semi-norms do not belong to L, because in 
these cases the constant C p may depend on the subdivision ao < ... < a n . These two 
examples are important in certain contexts, and this is why we would like to introduce a 
weaker definition, which will be given in terms of the evaluation of 1 1 • 1 1 along specific families 
of functions. 

If I £ N and tp £ C l K , we introduce \P = {ipj n , n £ j > 0}, the two-parametric subset of 
C l K defined by 

^j 3n {t)=^{t-n)) 

for every n £ Z, j > 0,t £ R. 

Definition 3 Let (3 € R ; p £ (0, +oo] and \ \ ■ \ \ be a lower semi- continuous, I -smooth- finite 
semi-norm on T . We say that \ \ ■ \ \ is a lower (/?,p)-semi-norm in the wide sense with respect 
to , if it satisfies the following property: 

(G) There exists C p > such that for every j > 0, arbitrary x±,..., x m £ M ; and arbitrary 
ni, . . . ,n m £ N such that the supports of the functions {ipj,^, 1 < * < m} have disjoint 
interiors, 

l/p 

if p < +00, 

,i=i J 



i=l 



< C p 2^ 



x i" l Pjn i 



i=l 



[0,1] 



[0,1] 




< 2^ sup \ Xi \ 
Ki<m 



if p = +00. 



Similarly we will denote by L(/3,p) the set of lower (f3, p)-semi-norms in the wide sense, 
and set L for the union of the L(/3,p)'s. Analogously, we define N(/?,p) = U(/3,p) PlL(/3,p) 
and set N for the union of the N(/?,p)'s. 



Remarks (a) In the sequel, will be either a family of sufficiently smooth wavelet func- 
tions, or some Schauder system on [0, 1]. 

(b) Condition (G) just means that (G) holds for linear combinations of specific functions 
with the same "frequency", and it is a well-known condition to obtain lower bounds for 
small deviation probabilities in a Gaussian framework [HHl [101 • Actually we will use it for 
the same purposes, but sometimes in a more general context, see Section @J 

(c) At first sight, neither translation-invariance nor self-similarity are required in the defi- 
nition of Jj(P,p). Actually, working with a specific family of function allows us to combine 
these necessary properties in the single inequality given by (G). 
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(d) Again we can prove that L(/3,p) C L(f3,p') (resp. N(/3,p) nN(/?',p') = 0) whenever 
p' < p (resp. ((3,p) / {(3',p')). 

(e) It is a bit tedious but not difficult to see that (rj,p,q)-Besov semi- norm belongs to 
N(ry, +00) and that (ry,p)-Sobolev semi-norm belongs to N(?7 — l/p,p), both with respect 
to the Schauder system. 

Notation In the following we will mainly consider functions restricted to [0, 1]. By minor 
abuse of notations, and for the sake of brevity, we will set ||/|| = ||/||[o,i] for every / G J 7 , 
although such statements like "|| • || G U(/3,p)" or "|| • || G L(/3,p)" will always refer to 
the family ||-|| = {|| - ||j, By C we will always mean a positive finite constant 

independent of the involved parameters, and whose value may change from line to line. 



2.2 Riemann-Liouville processes and their associated linear fractional sta- 
ble motions 

Let {Z t , t > 0} be a symmetric a-stable process with index a G (0,2], i.e. Z is a Levy 
process whose characteristic function is given by 



E 



e iAZ t 



e -t\\\° 



for every t > and A G M.. It is well-known and easy to see that for every H > the 
following stochastic integral: 



Rf = j\t-s) H - l l a dZ s , 
Jo 



is well-defined for every t > 0. We set R$ = and call {i?^, t > 0} the Riemann-Liouville 
process (in abridged form: RLP) with Hurst parameter H. This latter terminology is 
motivated by the following if-self-similarity property of R H : for every c > 

{R*, t > 0} i {c H R», t > 0} . 

When no confusion is possible, we will drop the subscript H and write R = R H for the sake 
of brevity. Notice that R has no stationary increments, unless H = 1/a. However R shares 
some kind of extrapolation-homogeneity, which will be important in the sequel: namely, if 
we set t 

Ra,t = [ (t-s) H - l / a dZ s 
J a 

for every t > a > 0, then the following equality in law holds: 

{Ra,a+t, t >0} = {R t , t>0}. 

Besides, R shares an equally important independence property: if R a j = Rt — Ra,t, then 
for every a > the processes {R atCl +t, t > 0} and {R a , a +t, t > 0} are independent. These 
three properties are easy to see as respective direct consequences (in reverse order) of the 
independence, stationarity, and stability of the increments of Z. 

Of course, Z itself is an RLP with Hurst parameter H = 1/a. Notice also that up to 
normalization constants, the family of Riemann-Liouville processes is closed with respect to 
time-integration. In particular, the m-times integrated Brownian motion is an RLP with 
parameters a = 2 and H = m + 1/2. 
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The Riemann-Liouville process is closely related to {X^ , t > 0}, the so-called linear stable 
fractional motion (in abridged form: LFSM) with Hurst parameter H. X H can be defined 
through the following (independent) decomposition: X H = R + M H where M$ = and 

Mf = £°° ((i + s) H - l / a - s H - l / a ) dZ s 

for every t > 0, Z being an independent copy of Z. We call M H the long memory process 
(in abridged form: LMP) associated to X H . Notice that the stochastic integral defining 
M H diver ges a.s. at +00 as soon as H > 1. Hence, X can be defined only for H € (0, 1). 
Again, when no confusion is possible, we will drop the subscript H and write X = X H 
(resp. M = M H ) for the sake of brevity. 

X is also ii-self-similar but, contrary to R, does not share the above extrapolation- 
homogeneity property. Instead, its increments are truly stationary: 

{X s+t -X s , t >0} = {X t , t>0} 

for every s > 0, and X is a so-called H-sssi process. We refer to the monographs |35| |17| 
for an extensive account on these latter processes. Notice that in the Gaussian case a = 2, 
X coincides with the canonical fractional Brownian motion, up to some numerical factor 
(see Proposition 7.2.6 in |35j). Moreover, an alternative definition of X can then be given 
through the following ("well-balanced") decomposition: 

x t = r°° (it+ S |*-v2_| S |H-iM 
j —00 

where B is the usual white noise (see Exercise 7.2 in 35 ). When a < 2, the corresponding 
integral 

/+00 . . 
(\t + s\ H - 1 / a -\s\ H - 1 / a ) Z{ds) (3) 

also makes sense, where Z is the so-called symmetric a-stable noise. However, X is no more 
equivalent to X (see Theorem 7.4.5 in [HUD- Actually, many other "bilateral" -ff-sssi a-stable 
processes can be constructed, which are all non-equivalent except in the Gaussian case (see 
Definition 7.4.1 and Theorem 7.4.5 in [HE])- Hence, our unilateral LFSM process is just one 
possible stable generalization of fractional Brownian motion. However, we will restrict our 
study to this unilateral process, even though our results can probably be adapted to some 
other situations. 

Notice that the process M is smooth on (0, +00), so that X and R exhibit similar local 
properties. In particular it is well-known (see Chapters 9-12 in |35| ) that R (resp. X) admits 
a continuous version if and only if 

a = 2 or H > l/a. 

In Section 6 we will give some other local properties of R and X, related to the examples 
of semi-norms (a)-(f) listed above. We finally refer to [55] for a more thorough comparison 
between R and X, in the context of econometric applications. 

3 Existence of the small deviation constant for RLP 

We can now state the main result of this paper. 



8 



Theorem 4 Let \\ ■ \ \ £ ~U(f3,p) and R be an RLP with Hurst parameter H. Assume that 
H > (3 + 1/p and set 7 = (H — [5 — l/p)~ l . Then there exists K. € (0, +00] such that 

lime 7 log P[||i2|| < e] = -K. (4) 

ej,0 

Remarks (a) It is interesting to note that the stability index a does not directly show up 
in the expression of the small deviation rate. 

(b) This result says nothing about the finiteness of the constant IC, which is of course a 
very important feature. To this aim, one needs to show a lower bound for small deviation 
probabilities with appropriate order - see Section 0] While a priori nothing indicates that 
the order 7 of small deviation probabilities we propose in Theorem |I] is the right one, it is 
indeed correct for all available examples, as soon as p is chosen as small as possible - see 
Section H3 

The next two paragraphs are devoted to the proof of Theorem |1] and we fix once and for 
all an upper (/3,p)-semi-norm || • ||, and an RLP R with Hurst parameter H > (3 + 1/p. 
We begin in proving a kind of stochastic superadditivity property, which will be our crucial 
argument. 



3.1 Stochastic superadditivity 

Let X and Y be two real-valued random variables. We say that X is stochastically larger 
than Y if and only if 

P{X < r} < ¥{Y < r} 

for every r € R, and we write in this case X >z Y. If X and Y are positive, then clearly 
X >z Y entails that 

E[exp-AX] < Efexp-AF] 

for every A > 0. 

Proposition 5 Suppose that p < 00. Let R\, R2 be two independent copies of R and set 
q = p(H — 0) > 0. Then for every a,b > 

{a + b) q \\R\\ P h a q \\R l \\ p + h q \\R 2 \\ P . 

Proof. Fix a, b > and set c = a + b. We first appeal to the p-superadditivity of 1 1 ■ 1 1 and 
get 

11-^1 1 [0,c] - 1 1^1 1 [0,o] + ll^ll[o,c]- 
By self-similarity for both || • || and R, we also have 

ll^llfcc] = H^aCc-Ml^ = ^^ll^llfo,!] = 

and, similarly, 

\\R\\\ , a] = a^R^ and ||22||J, j6] = b q \\R 2 \\P. 
Putting everything together, we now see that it is enough to show 

iWn 1 + ll^llf 1 >~ \\ R \\fn 1 + ll^allf 1, (5) 
11 n [0,a] 11 "[a,c] — 11 1 1 [0, (3.J 11 u|l [a,c]' \ ' 

because of the translation- invariance of 1 1 • 1 1 and the extrapolation-homogeneity of R. Here 
and throughout this section, R a denotes the process {R a ,t, t > a} For every a > 0, set Q a 
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for the filtration generated by {Z s , < s < a} and P a for the conditional probability with 
respect to Q a . We clearly have 



\R\\m i + 1 1^1 If i < r 
I N[o, a ] N M[a,c] — 



E 
E 



|i?||f , <r-\\R\\l , 

I 1 1 [a,cj — II n [0,a\ 



\Rn + Rn I If i < T — I \R\ |r n i 

1 " " 1 1 [a,cj — M M [Q,aJ 



for every r € R. Since i? a is a conditionally Gaussian process under P a (see p. 153-154), 
we can apply Anderson's inequality under P a and get (recall that R a is ^-measurable) 



\R a + Rn I If 1 ^ v — | |i?| If, , 

1 " u|l [a,c] — 11 "[0,a] 



< 



I D IIP ^ r _ I I T>\ IP 

l- fta ll[a,c] — r H-^NlCa] 



Now, since R a is independent of Q a , we can average backwards and obtain 



Wn I + IW 1 <r 

I <<[0,a\ n "[a,c] — 



for every r G R, as desired for 



< 



Wn 1 + ll^allf 1 < r 
I <<[0,a\ ii u n|a,c] — 



□ 



The next proposition covers the simpler situation when p = oo, and follows mainly the 
outline of Theorem 2.1 in |22j . 

Proposition 6 Suppose that p = +oo and set q = H — (3 > 0. For every a, b > and r G R 

P[(a + 6) 9 ||i2|| < r] < F[a"||-R|| < r] P[b 5 ||i2|| < r] . 

Proof. Fix a, 6 > and set c = a + 6. Again we appeal to the oo-superadditivity of 1 1 • 1 1 and 
get 

Il-R||[0,c] > SU P { I \R\ I [0,a], I \R\ I [a,c]} ■ 

Therefore, 



|i2||[o, c] < r} < P [||i?|| [0 , a] < r,\\R\\ M < r} = E 



l {\\R\\lO,a]<ry 



\\\\R\L,c] <H 



for every r € R, where P a is defined as above. By Anderson's inequality, independence of 
R a with respect to Q a , extrapolation-homogeneity of R and translation- invariance of [| • ||, 
we obtain 



1^1 l[a,c] < r} 



\Ra +-Ro||[a,cl < r 



< F a [\\Ra\\ M <r] = F[\\R\ 



[0,6] 



< r . 



Averaging backwards, we get 



P [||i% c] < r] < P [\\R\\ M < r] P [\\R\\ m < r] 
for every r G R. Now this yields the desired inequality 

P^ll^ll < r] < P[a 9 ||i2|| < r] W[b q \\R\\ < r] , 
since by self-similarity of R and 1 1 • 1 1 

P [||-R||[o,t] <r] = W[t q \\R\\ < r] 

for every t > and r G R. 



□ 



We are now ready to proceed to the proof of Theorem 0J 
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3.2 Proof of Theorem H 

Again we first consider the case p < oo. Proposition|5]yields the following decisive inequality 
for Laplace transforms: 

E[exp-(a + 6) 9 ||i2|| p ] < E [exp -a 9 p|| p ] E [exp -b q \\R\\ p ] 

for every a, b > and q = p{H—0) as above. This entails that the log-Laplace transformation 
of | \R\ | defined by 

*(/t) =logE[exp-/i 9 ||.R|| p ] 

for every h > 0, is a continuous negative function which satisfies ^f(a + b) < $?(a) + ^f(b) for 
every a, b > 0. By the standard subadditivity argument, we obtain 

]im m = M m = - c e [-00,0) 

h^cx) Ai /i>o n 
and, returning to the Laplace transform, 

lim A 1 /" E[exp-A||i?|n = -C. 

A— >+oo 

Notice that q > 1 by assumption. Hence we can appeal to de Bruijn's exponential Taube- 
rian theorem (see Theorem 4.12.9 in p], or Theorem 3.5 in [2^j for a more comfortable 
formulation), which yields 

lime 1 ^- 1 ) P[\\R\\ P <e] = -K = -{q - l)(C/q) q/{q - 1] , 

and finally 

lime 7 V[\\R\\ < e] = —K, E [— oo, 0) 

with 7 = (H — P — . This completes the proof of Theorem |I] when p < +oo. 

In the case p = +oo we do not even need Laplace transformation, since it follows directly 
from Proposition that the function 



q>(h) = -logP 
is subadditive. Again this entails 



\R\\<hP- H 



hm — ^ = mf — )-?- = —K E -oo,0), 
/i-too n, /i>o n 

and we obtain 

lime 7 P[||i?|| < e] = —JC € f-oo,0) 
with 7 = (i? — /3) _1 , as desired when p = +oo. 



□ 
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4 Lower bounds: finiteness of the constant for continuous 
RLP 

In this section we obtain a suitable lower bound for small deviation probabilities which will 
allow us to prove, under certain conditions, that the constant /C from Theorem 0] is actually 
finite whenever || • || G L(/3,p) as well. Unfortunately, our method is only efficient in the 
continuous case, i.e. when 

a = 2 or H > 1/a. 

An explanation for this understandable, but important restriction on R will be given later. 
Our result reads as follows: 

Theorem 7 Let || • || G L(/3,p) and R be a continuous a-stable RLP with Hurst parameter 
H. Suppose that H > (3 + 1/p if a = 2 and H > (3 + 1/p + 1/ 'a if a < 2. Then 

liminfe 7 logP[||i?|| < e] > -oo, (6) 

with 7 = (H-/3- 1/p)" 1 . 

Combining Theorem 0] and Theorem yields readily the following fairly general small 
deviation theorem for continuous Riemann-Liouville processes: 

Corollary 8 Let || • || G N(/5,p) and R be a continuous a-stable RLP with Hurst parameter 
H . Suppose that H > (3 + 1/p if a = 2 and H > (3 + 1/p + 1/a if a < 2. Then there exists 
K, G (0, oo) such that 

lime^logPOli?!! < e] = -K, 

e|0 

with 7 = (H - (3 -1/p)- 1 . 

We stress that if the parameter H is not too large, more precisely if H < 2, then it is pos- 
sible to obtain the lower bound for small deviation probabilities under a weaker assumption 
than || • || G L(/?,p). Here we just need that || • || G L(/3, p) with respect to 

f = {^ in = V(2 j i - n + 1), 1 < n < V, j > 0} , 

the Schauder system generated by the triangular function ip(t) = l[ 0j i](i)(l — |2i — 1|). 

Notice that on each level j the supports of the functions ipj n have disjoint interiors, so 
that here || ■ || G L(/3,p) w.r.t. ^ simply means that 



2 j 








n=l 





1/p 



for p < +oo, and with an obvious modification for p = oo. We have an analogous result to 
Theorem |7| 

Theorem 9 Let || • || G L(/3,p) w.r.t. Schauder system R be a continuous a-stable RLP 
with Hurst parameter H < 2. Suppose that H > (3 + 1/p if a = 2 and H>(3 + l/p + l/a if 
a < 2. Then 

liminf e 7 logP [||i?|| < e] > -oo, 

with 7 = (H - (3 -1/p)- 1 . 
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And, of course, we get the corresponding corollary: 

Corollary 10 Let \\ ■ || £ N(/3,j>) w.r.t. Schauder system and R be a continuous a-stable 
RLP with Hurst parameter H < 2. Suppose that H > /3+1/p if a = 2 and H > /3+1/p+l/a 
if a < 2. Then there exists K, € (0,oo) such that 

lime 7 log P[||i2|| < e] = -/C, 

with 7 = (H-/3- l/p)~ l . 



4.1 Decorrelating stable arrays 



In this paragraph we prove a crucial lower bound, which is a generalization of Lemma 2.1 
in |4U| to arrays of symmetric a-stable random variables, and which will be useful for both 
Theorems and El It relies on a version of Sidak's inequality for stable variables recently 
obtained in this framework by G. Samorodnitsky |34| . 

Lemma 11 Let M, h > and {yj n , 1 < fi < M2" J , j > 0} be an array of identically dis- 
tributed symmetric a-stable random variables. Let z, 5 > be such that 5 < z if a = 2 and 
5 < z — h/a if a < 2. Let m > be an integer and set 

dj = dj{m) = 2 z ^' m ^~ m \. 

Then there exists a constant C depending only on M, h, z, 5 such that 

\Vjn\ <dj, l<n< M2 hj , j > > exp -C2 hm . 

Proof. Up to normalization, the case a = 2 is just the statement of Lemma 2.1 in |4f)j . 
Hence we can concentrate on the case a < 2, and first notice that Lemma 2.1 in [51] entails 
the following decorrelation inequality: 

iM2 h i 



WjT 



<dj, l<n< M2 hj , j > 



ilL 

i>o 



y\ < dj] 



where y is some symmetric a-stable random variable. We can decompose the right-hand 
side into 



IL 



\y\ < dj 



M2 h i 




y\ < dj 



M2 h i 



m—1 



W^[\y\<d 3 ] 

j=0 



M2 h i 



To estimate the infinite product, we use the following well-known tail behavior of y (see e.g. 
Property 1.2.15 in 



lim r c 

rt+oo 



\ y \ > r] = K\ G (0,+oo). 



(7) 



This yields 



log JJ P [|2/| <dj 

j>m 



\M2 h J 



> 



> 



> 



M ^2^1og(l 

j>m 

-C 2hj 



D [|y|>^-]) 

y\ > 2^ _5 ^- 7_m ' ) 



j>m 

q c^j 2-oi(z~b){j-m) 



j>m 



-C2 



hm 
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where in the last inequality we used h — a(z — 5) < 0. The estimate of the finite product is 
even simpler. Since y has a positive density in the neighbourhood of the origin, we have 

lime _1 P[|j/| < e] = K 2 € (0,+oo). 

ej.0 



This entails 



m— 1 

log \{ T[\y\ <dj 

j=0 



m—1 



M 2 hj 'logP[|y| < dj] 



> 



3=0 
m—1 

C^2^(l + log(2( z+5 )( m ^)) 

3=0 

m—1 

> -C2 hm ^ (m - j) 2 h ( j - m 1 

j=o 

> -C2 



hrn 



where in the last inequality we used h > 0. Putting everything together now yields 

\Vjn\ <dj, l<n< M2 hj , j > oj > exp -C2 hm 
for a constant C not depending on m. 



□ 



Remark In the Non-Gaussian case, it is easy to see that the condition 5 < z — h/a is also 
necessary, because of the heavy tails of a-stable random variables: indeed, if the yjn's are 
mutually independent and if z = 6 + h/a, then it follows from (0) that 



lUjnl 



<dj, l<n< M2 hj , j > 



0. 



4.2 Some elements of wavelet theory 

The proof of Theorem [7| relies mainly on a suitable wavelet decomposition of R, which we 
recall here for the sake of completeness. In this paragraph we fix once and for all a semi-norm 
|.| € Jj(P,p) which is ^-smooth-finite, and R a continuous RLP satisfying the assumptions 
of Theorem [7| 

There exist two functions ip, ifi 6 C e K ("wavelet parents") such that tp has vanishing mo- 
ments up to order t: 

/oo 
t k ip(t) dt = < k < £, (8) 
-oo 

and such that the wavelet functions {ipj n , n G Z, j > 0} and {tp n , n € Z}, respectively de- 
fined by 

ip jn (t) = 2 j/2 ip(2 j t-n) and <p n (t) = tp(t-n), 

form an orthogonal base of L2(M). We refer to Daubechies' construction Section 6.4 
for the definition and Section 7.1 for smoothness properties of these compactly supported 
wavelets. A useful book is also [3U] . 
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Fix [-D, D] an interval containing the support of tp. Consider /, the integration operator 
on compactly supported functions: 



If(t) 



f(s)ds 



for every such / and t G M. We set 1° for the identity operator and I k for the fe-th iteration 
of J, k > 1. Since for every k > 



f(s)(t-s) k ds, 



the moment condition © on ift entails that the functions I k ip are also supported by the 
interval [— D, D] for < k < £. In particular, 



I k i(>(±D) = 0, < k < £. 



(9) 



It is well-known that under our assumptions on H, the process 1[ 0j _d+i]-R belongs a.s. to 
L2(M) (see |H5] Chapter 11). Hence, we can write its wavelet decomposition: 



E 

n 



I'n'-Pri i 



(10) 



where 

Tj n = / R(s)ip jn (s) ds and r n = / i?(s)(^ n (s)ds. 
io io 

Actually more can be said about the convergence of the series on the right-hand side in 

(|10|) . Namely, again from our assumption on H, it is well-known that R is locally r/-H61der 

for some r\ > |19| |41j . and in particular the function R coincides on [-D, 2D + 1] with an 

L2-function which is globally //-Holder. We know then (see e.g. Theorem 7 in |30j . Chapter 

6) that the series ()1()|) converges to R uniformly on [0, 1]. 

Besides, in (fTT)|) we can delete each ipj n and ip n whose support does not overlap with [0, 1], 

and the remaining series still converges to R uniformly on [0, 1]. Since ip and ip have compact 

support, this leads to a decomposition of R into a functional array of exponential size: there 

exists a constant M depending only on D such that 



j>l \\n\<M23 J \n\<M 

4.3 Proof of Theorem [7| 

We may (and will) suppose that £ > H. For every n £ Z we will set ipo n = 4> n and r§ n = r n 
for the sake of concision. By lower semi-continuity and triangle's inequality, we clearly have 



j>0 



r jn' l Pjn 

\n\<M2i 



Besides, since the parent functions have support in [— D, D], for each j > we can split the 
family 

{^n, |n| < M2?) 
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into 2D subfamilies such that in each subfamily the functions have supports with disjoint 
interiors. Suppose first that p < +00. Condition (G) yields 

ii^ii < ce( E M p feir) • (ii) 

j>0 \\n\<M2i J 

On the one hand, by homogeneity, ^-self-similarity, contr activity and smooth- finiteness of 
|| • ||, it is easy to see that for every j, ra, 

IhMI < C2<W)i (12) 

for some constant C independent of j, n. On the other hand, if j > 1, 

r jn = 2 j/2 I R{s)^{2 j s - n) ds 
-00 







= 2~ j/2 / R(2~ j s)tp{s - n) ds 
Jo 

1 r+00 
= 2-(" +1 /2)j / R( 8 )i/>(s-n)da 



= 2- { - H+1 Wi f R{s + n)4>(s)ds, (13) 

where in the penultimate inequality we appealed to the //-self-similarity of R. Plugging 
(TT2ll and JIEJ) into ((TTJ) yields 

||A|| < Cj2^ iH ^~ 1/p)j ( ^p \r' jn \], (14) 



j>0 



\n\<MV 



where r'- n stand for the renormalized wavelet coefficients of R, viz. 

= ^ r ;,i U. + „Ms )is itj >i. <i5) 

Jr 

Notice that one can write the wavelet coefficients from (|15|) in the following integral form: 
if j > 1, then r' jn = T n y jn , where 

r° = J (^J (u + n - s)+ 'ip(u) du^j ds 

H' = H — l/a, and {yj n , \n I < M2 J , j > 0} is an array of identically distributed symmetric 
a-stable random variables. We first aim to prove 

supr n < +00. (16) 

n 

To get this uniform bound, we first recall that ip has its support in [— D,D], so that the 
integral defining r" can be rewritten as follows: 

n+D / rD 







(u + n — s)^ ip{u) dii^j ds. 



16 



We cut the domain of integration over s into [0, n — 2D] and [n — 2D,n + D]. The first 
integral is given by 

l\ = J (j (u + n — s) tp(u) du) ds. 

We first transform 

(u + n — s) ip(u) du 

-D 

through £ successive integrations by parts. Recalling @, we see that each time the border 
terms / m ^(±.D), < m < £, vanish. In the end, we obtain 

i-a-2D / i-D \ a 
il = C J (j (u + n- s) H '- £ I e ip(u) du) ds 

m-2D 



pn—zu 

< C {D + n-s) c 

Jo 



+00 



< C / s a{ - H - e) - 1 ds, (17) 



D 



where in the first inequality we used I > H > H' . The integral on the right-hand side is 
clearly finite and independent of n. 
The integral over the second domain is given, after a change of variable, by 

r3D / r-D \ a r3D / r3D \ a 

ll = J [J (u + 2D- a)jV(u) duj ds < C J (J u 11 ' du) ds. (18) 

Since H' > —1/2, the integral on the right-hand side is again finite and independent of n. 
Putting (|T7|) and (|TH|) together yields (|Pd|) as desired. 
The following upper bound on ||i?|| is a direct consequence of f|14|) and (|TT))) : 

ll^ll < cVrW- 1 /* f sup | %n |), (19) 

j> \l<\n\<M2J J 

where {yj n , \n\ < M2 J , j > 0} is an array of identically distributed symmetric a-stable 
random variables on R. 

The end of the proof is now standard and follows 0Q], Theorem 3.1. Choose 5 > such 
that 

( 5 < H - (3 -1/p if a = 2, 

\ 5 < H - p -1/p-l/a if a < 2. 

Let m be a positive integer. Set 

dj = 2( // -' 3 - 1 /p)(i-™)-<5|j-m| (20) 
for every j > 0. On the one hand, it is clear from (|19|) and (|2U|) that 

{\y jn \ < dj, 1 < |n| < M2 j , j > 0} C < C(<5)2-^ /3 - 1 / p ) m } . 

On the other hand, it follows from Lemma fTTI (with h = 1) that 

P [\y jn \ <dj, 1 < \n\ < M2 j , j > 0] > exp-C2 m . 
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Since H — j3 — 1/p = I/7, we obtain 



liminf 2 _m logP [||JJ|| < C(5)2~ m ^ 



> —00 



which is equivalent to 



liminfe 7 logP[||i?|| < e] > -00, 

e|0 



as desired. The proof is complete in the case p < 00. The case p 
exactly in the same way, replacing by 



00 can be handled 



< cy;( 

j>0 \ 



sup \fjn\ 1 1 Yjn\ 
\n\<M2i 



□ 



Remarks (a) It is clear from the above proof that the assumption |.| € L(/3,j>) is not 
necessary to get the lower bound. We just need ||.|| G L(/J,p) w.r.t. ty, a wavelet family 
generated by parents {ip, ip} that are smooth enough. But the required smoothness depends 
on the parameter H, and this would lead to much heavier notations. For this reason we 
prefer stating Theorem [7| in this form, save for the loss of generality. 

(b) There is at least an example where L(/3,p) is really more relevant than L(/3,p). It is 
well-known (see [3U] Chapter 6) that (rj,p, g)-Besov semi-norms, rj > and p,q> 1, are 
equivalent to sequential norms on the wavelet coefficients. More precisely, 

n j,n 

Hence, with our notations, it is clear that the (r],p, g)-Besov semi-norm belongs to Jj(t] — 
l/p,p) w.r.t. any wavelet basis ^. Hence, when p = +00, Theorem 0] yields the existence 
of the constant for R with 7 = (H — rj) . Note that this rate is in accordance with the 
results of |4U| . which covered the range of parameters a = 2 and < H < 1 (and with no 
condition on p) . 

(c) In the Gaussian case, the smooth-finiteness of the semi-norm is obtained for free when 
||i?|| is a.s. finite. Indeed, if TL denotes the reproducing kernel Hilbert space associated with 
R, then it is well-known that Tt contains functions with [H + l/2)-th derivative in L2. On 
the other hand, from a 0-1 law for Gaussian measures on linear spaces, the finiteness of 
with positive probability yields that ||/|| < 00 for every / G TL. Hence, is ^-smooth- 
finite as soon as i > H + 1/2. In the non-Gaussian case, the challenging question whether 

< 00 a.s. implies the smooth-finiteness of the semi-norm remains open. 

(d) For the first efficient use of wavelet methods in similar problems to small deviation 
probability, we refer to where the optimal finite-dimensional approximation of fractional 
Brownian motion is considered. 
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4.4 Proof of Theorem EH 

The outline of the proof is the same as that for Theorem Q except that we use Schauder 
system and provide new estimates for the corresponding coefficients. Recall that from our 
assumptions, < H < 2 in the Gaussian case and 1/a < H < 2 in the non-Gaussian case. 
Again we set H' = H — 1/a. We exclude the case H' = where R is Brownian motion: the 
result follows then directly from Theorem 1 in [39J. 
Since R is continuous, we can decompose it on [0, 1] along the Schauder system: 

Rt = ^Hj{t) + Rit, 

for every f E [0,1]. Here, for every j > 0, 



2J 

£ 

71=1 



where we set t° jn = (n - 1)2^', t) n = (n - 1/2)2^', t 2 jn = n2~i and 



2R t i — R t o — R t 2 . 



We first suppose that p < oo. Condition (G) entails 



■v 



i/p 



\\HjW < cv* 2>i»n • (21) 

Vn=l 

Notice that the coefficients rj n can be rewritten as r, n = crj n yj n , where 



a 



J" 



/ +00 |2(4-,)f-(4-,)f-(4-,)f 



(22) 



(we use the notation m + = u A for every n G M), and {yjn, 1 < ra < 2 J , j > 0} is an array 
of identically distributed symmetric a-stable random variables. For each 1 < n < 2 3 , j > 0, 
we will now give an upper bound on <jj n depending only on j. 
If n = 1, then 



.A 



2 _J 



2 2 



-(j+l) 



(2-i - s) 



H' 



ds < C 

Jo 



aH ' ds = C2~ aHj . 



If n > 1, then we set i - = (n — 2)2 J > and cut the domain of integration in ()22j) into 
[0, tj n ] and [t~ n ,+oo). Reasoning as above, we see that the second integral 



.in 



2(t)n-sf' -(t° Jn -sf' -(t%- S f 



ds, 



is bounded from above by C 2 aH K We estimate the first integral as follows: 



1^ 2{t) n -sf -{t%-sf -(t%-sf 
Jo 



2~ j 



ds 



du 
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^' u H ' 



a 

du 



= 2~ aH i J 2{u- 1/2) H ' - (u- 1) J 

/■oo 

< C 2~ aHj / U ( H '-V a du 
= C 2- aHj , 

where the last equality comes from H < 2. Putting everything together yields now 



and, recalling l|21jl. 



CTjn < C2- Hj 



\TZjW < C2^ +1 ^ H ^ sup \y jn \ . (23) 

\ Kn<23 I 



By lower semi-continuity and triangle's inequality, this entails 

\\R\\ < C f \Ri\ + V2^ 1 / p -^ ( sup \y jr 

\ j>0 \l<\n\<M2i 

and since R\ has a symmetric a-stable law, we can finish the proof as in Theorem [7| 

□ 



Remarks (a) In the latter proof, the estimate on lj n becomes too crude when H > 2. 
This is the principal reason why we need to introduce smoother wavelets in Theorem 

(b) As we said before, the proofs of Theorems and El work only for continuous processes. 
The main reason for this comes from Lemma 1111 On the one hand, the exclusion of the 
boundary value z = 5 + h/a in the non-Gaussian case cancels important case H = l/a>. On 
the other hand, Sidak's inequality makes us handle the ||7£j||'s as if they were independent. 
But if R is discontinuous, then it is possible that its jumps have a significant influence on 
every level IZj, so that our estimate of the series Y2j ll^jll m &y not be realistic anymore. 

(c) In the Non-Gaussian case and when H < (3 + 1/p + 1/a, it is easy to see that the 
statements of Theorems [7] and 03 are false, even if H > (3 + 1/p. Suppose for example that 
H < 1/a. It is well-known (see [SSI Chapter 10) that a.s. 

sup \Rt\ = +oo, 

0<t<l 

so that Theorems and |H1 cannot hold with respect to the supremum norm (i.e. /C = +oo 
in Theorem|l]), although here (5 = and p = +oo. Similarly, if H = 1/a, then clearly 

(\Rt — Rs I \ 
— ; j=- = +oo a.s. 
\t — s\ J 

for every n > 0, so that Theorems and El cannot hold with respect to any ^-Holder semi- 
norm (f3 = rj,p = +oo). It is natural to conjecture that these two theorems remain true in 
full generality when H = /3 + 1/p + 1/a, but in view of the above Remark (b), we probably 
need different methods. 
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5 Lower bounds for LMP: Finiteness of the constant for con- 
tinuous LFSM 



Our aim in this section is to extend Theorem|§]to X, the continuous LFSM which we defined 
in Section 2. Recall that X admits an independent decomposition X = R + M where R 
(resp. M) is an RLP (resp. an LMP) with the same parameters. The following theorem 
can be viewed as a generalization of Lemma 3.2 in |22| . 

Theorem 12 Let || • || £ L(/?,p) w.r.t. the Schauder system and M be an LMP with 
parameters a £ (1>2], H £ (1/a, 1). Suppose that H > j3 + 1/p if a = 2 and H > 
p + 1/p + 1/a if a ^2. Then 



where 7 = (H — j3 — 1 /p) 1 . 

Using Theorems 0] and |H1 along with the elementary independence argument developed in 
22 p. 1334, Theorem 1121 yields readily the desired small deviation theorem for X: 

Theorem 13 Let \ \ ■ \ \ € N(/?,p) w.r.t. the Schauder system and X be a continuous LFSM 
with parameters a £ (1,2], H £ (1/a, 1). Suppose that H > f3 + 1/p if a = 2 and H > 
(3 + l/p + 1/a if a ^ 2. Then 

lime 7 logP[||X|| < e] = lime 7 log P < e] = -K 

sj.0 ej.0 

where JC € (— 00, 0), 7 = (H — (3 — 1/p)' 1 and R is the RLP associated to X. 

We now proceed to the proof of Theorem 1121 which is quite analogous to the proof of 
Theorem |UJ except that our estimates on the Schauder coefficients will not be uniform in 
time-argument since M has a singularity at 0. 

5 . 1 Proof of Theorem fT2l 

Again we set H' = H — 1/a £ (—1, 1) and exclude the trivial situation when H' = 0. Since 
M is continuous, we can decompose it on [0, 1] along the Schauder system: for every t £ [0, 1] 



lime 7 logP[||M|| < e] = 0, 



M t = ^Mj{t) +Mit 



where 




n=l 



and 

wijn — 2Af+i — ■^iP — -^t 2 ^ ^ th.6 proof of Theorem |^| We first consider the case 



I ■ L ■ 

jn jn jn 



p < +00. Condition (G) entails 




(24) 
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and we just need to evaluate the coefficients uij n — crj n yj n , where 



a 



j a 



2(s + t} n ) H '-(s + t° jn r -(s + t%) 



o \H' 



ds 



-Haj 



2(u + n- l/2) H ' -{u + n- l) H> -{u + n) H ' 



du, 



and {yjn, 1 < n < 2 J , j > 0} is an array of identically distributed symmetric a-stable ran- 
dom variables. Notice first that from the so-called Cooper's formula 



1/2 ru+n-l/2+8 
Ju+n-1+8 



2{u + n- l/2) H ' -(u + n- if - (u + nf = -H'(H' - 1) 

Clearly, since H' < 2, this entails 

2(u + n-l/2) H ' - (u + n-l) H> - (u + n) H ' < Cmin {(it + n - l) H '~ 2 , l} . 
We now fix some h E (0, 1) and first estimate aj n when 2 h i < n < 2 J : 



v H '~ 2 dvd6 



(7 



f + OD 

< C2- Haj {u + n-l) a(H '~V du 
Jo 

< C2- H ^(n-l)^ H '-^ +l 

Qt^—Hoiji^—hai^l—H^j 

The estimate of <jj n when 1 < n < 2 h i is even simpler: 



af n < C2~ Haj 



U 



°<H'-2) du+ / Cdu ) < C2- Ho *. 



(25) 



(26) 



Recalling we get 



/ 2« V \ 1/P 

\\MjW < C2« E°in+ £ .-p .V, 

\n=l n=2^'+l / 



l<n<2J 

;\1/P 



< C2 -(^Y 2 ^ + 2J2 -M2^ te y p sup |y . r 



The balance is attained at h = (1 + (2 — -£f)p) 1 £ (0) 1)) whence 

< C2~^ H - p - h ^ j sup |y in |. 

l<n<2^ 

Using (|27|) instead of (|23j) , we can proceed as in the proof of Theorem |U] and obtain 

liminfe 7 'logP[||M|| < e] > — oo, 

e|0 

with 7' = (i? — (3 — h/p)~ l < 7. In particular 

liminfe 7 logP[||M|| < e] = 0, 

ej.0 



(27) 



which completes the proof of Theorem ^] for p < 



00. 
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The case p = +00 requires more careful estimates. Condition (G) entails 

\\Mj\\ < C2 Pj sup \m jn \. (28) 

l<n<23 

Again we fix some h € (0, 1). Using (|28[). (]25[). and (|2*H|) we can write, 
||A4j|| < C2 ft sup Oj- n |y 3 - n | 

l<n<23 

< C , 2 -(H-/9)i[ gup j y . n | +2 -M gup 

\l<n<2>v 2 h i<n<2i J 

where we set hi = (2 — -£f)/i. For every integer m > we focus on the event 

= {\Vjn\ < dj n (m), 1 < 71 < 2 j , j > 0} 

where we set mi = [m/h], 5 E (0, H — (3) and 



djn(m) 



2(- H '^/ 3 )0'^ m i)-< 5 li- m il if 1 < n < 2 h i 

2(H-/3+h 1 )(j-m)-S\j~m\ •£ 2 h 3 < n < 2 J . 



Take now h small enough such that 

H - (3 < H - 13 + h x < (H - l3)/h. 
On the one hand, it is clear that 

n m C jx^ll-Mjll < C2-< H - /3+ ' ll > m j . (29) 

On the other hand, it follows from decorrelation argument of Lemma ITT1 that 

logP[O m ] > - C (2 hmi + 2 m ) > -C2 m . (30) 
Using (|29|) and (|3Ujl . we can now finish the proof exactly as in the case p < +oo. 

□ 

Remarks (a) Contrary to the above proof does not require any entropy argument 
and relies only on an elementary estimate of the Schauder coefficients. 

(b) It would be quite interesting to calculate the optimal rate 

7o = infi7>0 lime 7 logP[||M[| < e] = 

We did not study this question in details but it seems plausible that 70 = 0: indeed, M is 
a C°° process on (0,+oo), and one may expect sub exponential rate for its small deviation 
probabilities. 
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6 Examples: particular semi-norms 



In this section, we place all our results in the context of the previous literature, and show 
what is new. We try to be as exhaustive as possible as far as RLP's and LFSM's are 
concerned. We refer to the surveys j^S] |26| for more information about other processes and 
for further references. 

Everywhere R (resp. X) will be an RLP (resp. an LFSM) with Hurst parameter H > 
and stability index a S (0,2], whereas |.| will be a semi-norm in N(/3,p) or N(/3,p) for 
some (5 € M and p € (0, +oo]. We always suppose that H > (3 + 1/p and set 

j = (H-P-l/p)- 1 

for our small deviation rate. K, will stand for the small deviation constant appearing on 
the right-hand side of Theorem 01 Sometimes we will call them just "the rate" and "the 
constant" . Let us begin with the most classical case, which deserves of course a particular 
mentioning. 

6.1 Brownian motion 

Brownian motion B is an RLP (or an LFSM) with parameters H = 1/2 and a = 2, hence 
it clearly satisfies the assumptions of Corollary (or Theorem I1MJ1 . Notice that here our 
Theorem |§] just amounts to Stolz's lower bound criterion (Theorem 1 in [HH]); s ° that in the 
present paper the novelty comes only from Theorem 0J Corollary entails the existence 
of Brownian small deviation constants under almost all the classical semi-norms (with the 
regrettable exception of certain Besov semi-norms for which we were unable to prove the 
appropriate super additivity index). In other words, we get 

lime 7 logP[||5|| < e] = —K € (-oo,0), 

ej.0 

where the dependence of 7 on the semi-norm 1 1 • 1 1 is given by the following table: 





7 


Supremum 


2 


L p 


2 


7?-H61der 


2/(1-27?) 


p- variation 


2p/(p - 2) 


(7?, p)-Sobolev 


2/(1-27?) 


(77, 00, g)-Besov 


2/(1-27?) 



The existence of a finite small deviation constant for Brownian motion under L p and 
supremum norms is a very classical result, and in this situation the constant K, can even 
be more or less explicitly computed (see |2,~>| and the references therein). The existence of 
the constant under Holder semi- norms was first proved by Baldi and Roynette [2 ] , and it is 
actually a direct consequence of Stolz's upper bound criterion (Theorem 2 in [23) and of 
the standard subadditivity argument involving the Markov property of Brownian motion. 
In the three last examples, Stolz's results do provide the right rate, viz. 

-00 < liminf e 7 logP[||5|| < e] < limsupe 7 logP [||£|| < e] < 
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with the above corresponding 7, but here the existence of the limit does not follow directly. 
In these three cases, our result seems to be new. 

It is of course quite tantalizing to compute these new Brownian constants. As a rule, this 
kind of (hard) problems requires analytical techniques depending heavily on the choice of 
the semi-norm, and it seems difficult to tackle them through a global study. 

We stress that the above table matches quite accurately the sample path properties of 
Brownian motion. It is well-known that B is 77-Holder if and only if 77 < 1/2 and has finite 
p-variation if and only if p > 2. Besides, it has finite (r/,p)-Sobolev semi-norm if and only if 

77 < 1/2 na. 

Besov semi-norms deserve a special remark since no necessary and sufficient conditions 
seem to be available. Roynette jH2j had proved that Brownian motion has finite (f],p,q)- 
Besov semi- norm whenever one of the three following situations occurs: 



V 


p 


q 


< 1/2 


> i/v 


> 1 


< 1/2 


= 1/7? 


= 1 


= 1/2 


> 2 


= +00 



Notice that the latter situation is rather specific, since here small deviation probabilities 
vanish when e is small enough, thus the investigation is pointless. 

We are able to prove the existence of a finite small deviation constant in the following 
situation: 



V 


p 


q 


< 1/2 


= +00 


> 1 



6.2 Gaussian fractional processes 

In this paragraph a = 2 and H > 0, the case H = 1/2 being implicitly excluded. Up to 
normalization constant our process X is the usual fractional Brownian motion, while R can 
be viewed as the fractionally integrated Brownian motion. Before stating our results, we 
will give a list of local properties of R, most of which are classical. As we mentioned in 
Section 2, all these local properties remain true for X as soon as H < 1. 

R is //-Holder (resp. r?-Lipschitz) if and only if H > i]. In particular, R has finite p- variation 
if H > 1/p and only if H > 1/p (the boundary value being excluded when H > 1/2 - see 
Theorem 5.4 in [T^], but the situation H = 1/p < 1/2 seems to have escaped investigation). 
Finally, R has finite (r?,p)-Sobolev semi-norm if and only if H > 77 13 . 

Less is known about Besov semi-norms. In |12| . Proposition 4.1 and Theorem 4.3, it is 
proved that if H < 1 and < 1/p < rj, then X has finite (ry,p, oo)-Besov semi-norm if and 
only if H > rj (when H = rj, the small deviation probabilities vanish for small e). 

With the same notations as in the preceding paragraph, the following table is a direct 
consequence of Corollary |H1 and of Remark (b) at the end of paragraph 4.3. Recall that it is 
available for R and for X with the same finite constant IC, as soon as H > j3 + 1/p. 
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II ■ II 


7 


Supremum 


1/H 


L p 


1/H 


r/-H61der 


1/(H- V ) 


p- variation 


p/(Hp - 1) 


(r/, oo, g)-Besov 


1/(H- V ) 



If in addition H < 2, then we can add another line: 





7 


(?7,p)-Sobolev 


l/(H- V ) 



The existence of K, for X under the supremum norm was established in [22] and (3Z1- I n 
|22| . it is also proved for R when H < 1. When H > 1, the existence of /C is mentioned 
in [23] • Under Lp-norms, the right rate for i? can be found in j^j. Specialists seem to be 
aware of the existence of the constant j2H]) but up to our knowledge this is not explicitly 
proved in literature. Notice that the case p = 2 is rather special because of its Hilbert 
space background. Quite recently, the exact value of /C under the L2-norm was found for 
integrated Wiener process JOJ and for fractional Brownian motion |Sj. 

In all other situations, our results seem to be new. 
6.3 Non-Gaussian stable Levy processes 

Non-Gaussian symmetric a-stable Levy processes are RLP's with parameters a < 2 and 
H = 1/a. They have a.s. discontinuous sample paths, hence in the present paper they are 
only concerned by Theorem [I] 

Recall that they have finite p- variation if and only if p > a 6 j. In the case 1 < a < 2, 
it is shown in ^2] Theorem 6.1 that they have finite (1/a, p, oo)-Besov norm if and only if 
1 < p < a. 

Our Theorem |1] yields the following table which can be read as in the preceding paragraph, 
except that here we know nothing a priori about the finiteness of /C. 





7 


Supremum 


a 


Lp 


a 


p- variation 


apj (p — a) 



Under the supremum norm, the finiteness of /C dates back to Taylor [12] and Mogulskii 
|31| . Notice that in this case the existence of K, follows readily from the Markov property, so 
that our Theorem 0] is a bit useless here. Under L p -norms, the existence and the finiteness 
of IC seems to belong to some ancient folklore, but we could not find any precise reference in 
the literature. We remark that this follows readily from Theorem and Taylor- Mogulskii's 
lower bound. Z. Shi also indicated us that in this case the problem of L p -small deviations 
can be reduced to some specific large deviation probabilities for occupation measures which 
were previously studied by Donsker and Varadhan. Finally, it is shown in 38 that the case 
p = 2 is again rather specific since it can be derived from the L p -small deviations of Wiener 
process, thanks to an identity in law due to Donati-Martin, Song, and Yor |15j . 
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Nothing is known about the finiteness of /C for the p- variation norm, but there is no reason 
to believe that the rate 7 indicated above is not the right one. This will be the matter of 
further study. 

6.4 Non-Gaussian fractional processes 

In this paragraph a < 2 and H > 0, the case H = 1/a being implicitly excluded. As in 
the Gaussian framework, we begin with listing local properties of R, most of which can be 
found in the monograph [3S] . Again, all these local properties remain true for X when it is 
well-defined. 

R is sample-bounded iff H > 1/a and belongs locally to L p if and only if H > 1/a — 1/p. 
R is 77-Holder (resp. r^-Lipschitz) if and only if H > rj + 1/a. This latter fact is less classical 
and was proved for X by Takashima JHj (see also |19|). Notice that the "if part is quite 
straightforward from 

\Rt~Rs\ < C \\Z\loo |£ _ s |mm{l,.ff-l/a} 

for all s,t £ [0, 1]. It is interesting to note that contrary to the Gaussian case, the boundary 
value rj = H — 1/a is included here. In particular, this entails that R has finite p- variation 
if and only if H > 1/p + 1/a [II]. 

Nothing seems to be known about Sobolev or Besov semi-norms. From Holder continuity, 
we see that R has finite (r/,p)-Sobolev semi-norm if H > 77 + 1/a, yet this does not seem to 
be a sharp estimate. 

When R is continuous and when H > (3 + 1/p + 1/a, we get exactly the same tables as in 
paragraph 6.2. 

The only result which we are aware of in this direction is due to G. Samorodnitsky |34j . 
who had given just a bit less precise bounds in the small deviation problem for LFSM's 
(1/a < H < 1) under the supremum norm. In the present paper we find the right rate, 
prove that the constant exists, and our condition on H is the best possible since R is 
not bounded anymore when H < 1/a. In general, our results match quite accurately the 
local sample-path properties of R. Namely, except for Sobolev and Besov norms, the only 
situations which are not covered are the critical values of H w.r.t. Holder and p-variation 
semi-norms. 

When R is not continuous, i.e. when H < 1/a, the only relevant classical semi-norms are 
the Lp-norms, as soon as H > 1/a — 1/p. Our Theorem [1] yields 7 = 1/H but nothing is 
known about the finiteness of the constant JC. We believe that 1/H is still the right rate in 
this ultimate situation. 

7 Concluding remarks 

(a) In general, it is not true that H is the only parameter of a self-similar stable process 
which is involved in the exponential behavior of its small deviation probabilities. Indeed, 
|34j delivers an example of an a-stable, H-sss\ process whose small deviation rate depends 
both on H and a. Let us also recall the general results of «33] about stable measures, which 
entail that 

a 

7 < z 

1 — a 

for arbitrary a-stable processes with < a < 1. 



27 



(b) The symmetry assumption on the stable processes is a crucial one, since our proofs are 
heavily based on two Gaussian results: Anderson's and Sidak's inequalities. Nevertheless, 
it seems likely that everything remains true as soon as Z is not a subordinator. 

(c) Theorem |I] can be easily extended to multi-dimensional RLP's having some symmetry, 
since they also share the crucial extrapolation-homogeneity property. However, adaptating 
the proof of Theorem [7| in a multi-dimensional context raises some technical difficulties 
related to Sidak's inequality. 

(d) As we indicated before, there are many different generalizations of fractional Brownian 
motion in non-Gaussian case. If X is the well-balanced linear fractional stable motion 
defined in ©, then the proof of Theorem |I] fails in the absence of extrapolation- homogeneity. 
However, we believe that the rates should be the same. 

(e) Although our results are designed for nice translation-invariant semi-norms, it is worth- 
while to mention that there exist elaborated techniques concerning norms which do not share 
any translation- invariance. See especially |27| . |28| and |38| . where weighted L p -norms are 
handled respectively for Brownian motion, fractional Brownian motion, and symmetric te- 
stable processes. 
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